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We study the thermoelectric transport in two dimensional topological system which has co¬ 
existence of superconductivity(SC) and spin-density wave(SDW). The SC is presumed to be of 
d„, 2 _y 2 {px -\-ipy) type whereas the SDW order parameter is of BCS symmetry. The Hamiltonian 
describing such a coexistence phase is shown to have topological phase in addition to the conven¬ 
tional one. The transport properties in such topological system have two distinct contributions: 

(i) the surface/edge and (ii) the bulk. The competition between the surface/edge versus the bulk 
transport is analyzed in different parameter regimes and the possibility of enhancing the figure of 
merit is discussed. 
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In last few years, a new field in condensed matter sys¬ 
tems has emerged which is called ” Topological insulators 
and superconductors” Hi- It is based on the realization 
that the spin-orbit interaction in materials can lead to 
such electronic phases which has been observed in 
real materials^ topological insulator, similar to an 
ordinary insulator has a bulk energy gap separating the 
highest occupied electronic band from the lowest empty 
band [Tol - [T^ . It is closely related to the two-dimensional 
integer quantum Hall states [l3| which have unique edge 
states. The surface of such insulators however, neces¬ 
sarily has gapless states (edge states) that are protected 
by time reversal symmetry. These states are conducting 
with properties unlike any other known one-dimensional 
or two-dimensional electronic systems. These states are 
predicted to have special properties which are thought 
to be useful for applications ranging from spintronics 
to quantum computationfisjl . Due to such technologi¬ 
cal importance the field of topological systems has grown 
rapidly over the years. Since one of the compelling cri¬ 
teria of a material to be topological is that the system 
ought to have an energy gap, the topological concept ex¬ 
tends over to any gapped sterns such as superconduc¬ 
tors and superfluids etc. [IMI. Recently, a possibility of 
a topological phase in the coexistence phase of supercon¬ 
ductivity and spin-density wave has been discussed [T^ 
[2lj . The nature of such topological phase depends on the 
symmetry as well as the amplitudes of both the order pa¬ 
rameters. 


Recently, the photoemission and the scanning tunnel¬ 
ing spectroscopy have become an important experimen¬ 
tal tool in identifying the surface states in topological 
systems but the signature of such states in transport 
measurements is still under debate. The focus of ex¬ 
perimental studies thus has shifted to the separation of 
bulk conduction from the surface conduction. Recently, 
there is an attempt to formulate the thermoelectric trans¬ 
port in topologial insulators theoretically . Thermo¬ 
electric transport is basically conversion of heat to en¬ 


ergy. The efficiency of such energy conversion depends 
on the ’’figure of merit” of the material which is de¬ 
fined as, ZT = ^ where T is the temperature, 

S is the Seebeck coefficient and a and k are respectively 
the electrical and thermal conductivities. In order to 
achieve a high ZT, the system should be a good elec¬ 
tronic conductor whereas a bad lattice conductor. Fur¬ 
ther, the material low-dimensionality also have large S 
and have high efficiency due to their peaked density of 
statesjl^. Despite such proposals it was very hard to 
discover good thermoelectric materials in the past. This 
is due to the fact that a, S and k cannoot be indepen¬ 
dently controlled, namely, a material with large electri¬ 
cal conductivity a has large thermal conductivity n. The 
discovery of topological systems gives some hope on the 
above issues due to the fact that the edge state conduc¬ 
tion remains good whereas the phonon conduction is sup¬ 
pressed. The edge states are one-dimensional which satis¬ 
fies the low dimensionality [2^ criteria mentioned above. 
The topological phenomena which has already been ob¬ 
served in materials such as Bi2Se3^^ and 

Bi2Te3^^ are known as good thermoelectric materials. 

In this communication, we study the thermoelec¬ 
tric transport in two-dimensional topological system 
which has coexistence of superconductivity(SC) and spin- 
density wave(SDW). The symmetry of SC is considered 
to be of dx‘^_y2 -I- {px + ipy) type whereas the SDW order 
parameter is of Re'S" symmetry. The Hamiltonian having 
such a structure is shown to have topological coexistence 
phases in addition to the conventional one. The thermo¬ 
electric transport in such topological system is discussed 
with respect to two distinct contributions: (i) the surface 
and (ii) the bulk. The competition among the surface 
versus the bulk transport is analyzed in different param¬ 
eter regimes and the possibility of enhancement of figure 
of merit is discussed in these systems. 

In the present work, we consider the coexistence of 
SDW and d-wave superconductivity which can generate 
a triplet and non-zero center of mass superconducting 
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order parameter. We, thus start with a Hamiltonian on 

I 


a 2D square lattice [2l|, as, 




U 






r\ .ru\r\ 


V^{k, /c')4_.fcL^, 

k,k' 


k,k' 


( 1 ) 


Here, is the bare dispersion due to the tight 
binding approximation on a 2D square lattice, U 
is the on-site Coulomb interaction, are the 

pairing strengths for d-wave and p-wave super¬ 
conductivity and N is the number of sites. Also, 
^k = —2t{coskx + cosky) — it'cos kx cos ky — fi and 
4 (t (cfccr) denotes creation (annihilation) operator of 
the electron with spin a = (t,4-) at k = {kx,ky). 

Here, is the sum of k over the reduced Brillouin 
zone (RBZ). We express the wave-vector k in units 
of with ’a’ the lattice parameter of the underlying 
square lattice. Q = (tt, tt) is the SDW nesting vector 
in 2D. We assume here a commensurate SDW so that 
k + Q = k — Q. The staggered spin magnetization is 
defined as Mq = —^ J2k o- ^ ^ 4 +q a'^k,a >■ Since the 
discussion would be about three order parameters below, 
the crystal symmetry are such that the commutator of 
any two of them should give the third one. Thus, if 
is assumed to be of singlet d-wave symmetry, the 
SDW state guarantees that should be of triplet type. 
So we get the singlet interaction = VQSkSk' and 

= VgPkPk', where Sk = ^{coskx — cosky) and 
Pk = sin kx+i sin ky. We further assume that are at¬ 
tractive. The SC order parameters are defined as, for sin¬ 


glet state, Al, = A^sk' = Vg^Sk'J^k^k < >= 

< cl+Q,Y-k+Q,i >■ On the other hand, 
the triplet order is A^, = A^pk' = V^pk'Y^kPk < 
>= A2(sinA:,, -h zsinfcy) = Al^^ + iA\^ 
while All = V^pk>J2kPk < cl Y-k+Q,i >■ Defining 
= —4t' cos kx cos ky — p and = —2t(cos kx + cos ky) 

and employing the nesting property in the band disper¬ 
sion i.e. = C> C+Q = -C and also the order 

parameters A^_,_q = _= —Al and 
Afe+g = ~A|, the Hamiltonian in the momentum space 
can be expressed as, Ti = iljl'H(k)ipk where the four- 
component spinor V'fe is, V'l = (Cfct A-fc-Q4.> cLfcj.) Cfcs-Qt)- 
Thus, the Hamiltonian matrix 'H{k) in this basis is writ¬ 
ten as. 


t + ^k 

Ai 

Al 

Mq \ 

Af 


Mq 

-Al 

A^ 

Mq 

-(Cfe+4“) 

-Ar 

Mq 

-A^ 

-Al 
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In what follows, we study the energy spectrum of the 
above Hamiltonian. The Hamiltonian (eqn.(3)) is diago¬ 
nalized and the quasiparticle spectrum is obtained as. 


E±Mk) = + (Ai)2+ I Al P +M2 ± 2^4-^ I Al P +(Af^, Al -Mo^+)^ + 


-2 

k * 


It is obvious that the energy spectrum is fully gapped 
as shown in Fig. 1 and the gap closes only when the 
rhs of the above equation vanishes. A straight forward 
calculation provides a condition where the gap closes at 
points k = ( 0 , 0 ) and ( tt , tt ) is 


16t^ + = {it' + fi)^ . (3) 

Thus, one can find a topologically trivial and non-trivial 
regions with respect to these parameters. Since the crit¬ 
ical phase lines (for t = 0) are determined by the condi¬ 
tion Mq = {it' + p)"^ the phase becomes topological in 
the region where Mq < {it' + p)"^ whereas it is trivial for 


Mg > (4t' + p)"^- The details of these work have been 
discussed in a recent paper (2^. 

In order to study thermoelectric transport in these sys¬ 
tem, we consider a 2-D sample as a ribbon geometry 
where the ribbon width is taken to be very narrow. This 
is due to the fact that the edge states can have com¬ 
parable contribution as compared with the bulk. Using 
a linear response theoryj^, the electric current j and 
thermal current w which are coupled, are written as 


j/q 

w 


Lq Li 

Li L 2 



( 4 ) 


where q is the electron charge — e, p is the chemical po- 
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FIG. 1. Energy spectra E±^+{k), corresponding to coexis¬ 
tence of SC order parameters d^ 2 _y 2 -I- (px + ipy) and that 
of the SDW order parameter showing fully gapped spec¬ 
trum. For illustration, here, we have chosen t' = —0.3z, 
Mo = 0.0752, Aj = Ag = 0.0052, (2 = 0.30 eV). 


tential. Here, the field used are the electric field and the 
thermal gradient. The electric and thermal transport co¬ 
efficients are obtained as, 


a = e^Lo, S 


eT Lo’ 


zT = —r 


_3 _ 

Lf + klTLq 


1 L 0 L 2 — L\ 
T Lo 


(5) 


where Ke is the electron thermal conductivitiy and kl is 
phonon thermal conductivity. Here, L^’s are the corre¬ 
lation functions determing the thermoelectric transport 
coefficients. 

Since the topological systems have distinct edge and 
bulk states, we consider the transport due to both inde¬ 
pendently. We first consider the edge transport only. To 
describe the coherent transport of the edge states, we use 
the Landauer formula. The edge states are assumed to 
be perfectly conducting and the transmission coefficient 
T{E) is taken as unity. This is true when the electron 
energy is within the bulk gap (—A < if < A), A being 
the effective gap. The energy here is measured from the 
bottom of the conduction band. While carrying out the 
calculation of the transport coefficients, we consider the 
bottom of the bulk conduction band and neglect the va¬ 
lence band such that we deal only with the edge states. 
We further, restrict the chemical potential p to be within 
the gap. The edge state correlation faction L^, is given 
by 


where the suffix ’e’ means the edge transport, h is the 
Planck constant and £ and s respectively are the length 
and the cross-section of the sample, the factor 2 comes 
from the two gapless channels of the 2D system. The 
integral in the above correlation function can be rewritten 
in a dimensionless form as 


Ll{il) = -{kBTr nr 




(7) 


where y = ^, A = ^, /I = ^ and /3 = ^, T 
being the temperature. Using this edge-state correlation 
functions the thermoelectric coefficients are calculated. 
It is noticed that the edge contribution to the thermo¬ 
electric coefficients such as the figure of merit (ZT) is 
independent of the system size i and s eventhough the 
edge correlation functions are directly proportional to £ 
and inversely proprtional to the cross-sectional area s. 
Further, ZT is unusually large and exceeds unity when 
the chemical potential is in the bulk band. This is due 
to the fact that when y, is in the bulk band, the correla¬ 
tion functions Lq has only exponential dependence on y 
but Li and L 2 have combined exponential and algebraic 
dependence. On the contrary, when the chemical poten¬ 
tial is in the bulk gap, ZT becomes exponentially small 
{y -P 0) which is in agreement with the figures (see Fig. 
2 and Fig. 3). 

Next, we consider the bulk thermoelectric transport. 
The bulk transport are calculated using the Boltzmann 
transport theory. We assume that the relaxation time Te 
is constant. Since these transports are valid within the 
inelastic scattering length, we regard the inelastic scat¬ 
tering length as the effective system size. We use only the 
upper two subbands because there is a large gap between 
the upper and lower subbands due to the narrow-ribbon 
confinement. The bulk correlation function Ljy is written 
as 


Lliy) = J dEiE - yf D{E)Tev\ (8) 

where the suffix ’6’ means the bulk transport and Te is the 
relaxation time which is assumed to be constant. Also 
D[E) and v respectively are the bulk DOS and the veloc¬ 
ity. The above integral can be changed to integral over 
the momentum variables as 


lI{l) = 


cPk fdE++{k) 


(- 


J (27r)2 \ hdkx 

df{E++{k)-y) 


{E++{k)-yr 


dE^ 


) + (f;++^a+_). (9) 


-s-i- 


Using dimensionless variables namely kxO, = yi and 
kya = 2 / 2 , the above equation is written as 


= ^(kBTr+^£ dy, £ dy2{^^£^r 

p{L+{vi,V2)-p.) 

{E+{yi,y 2 ) - y) ^^^^E+(yi,y 2 )-p.) 

The variable in the above equations is written as, 
E+ — FiT which is functions of yi = kxa and ?/2 = kyO- 
Using the parameters in the energy dispersion, the above 
correlation functions are calculated. It is noted that all 
the bulk correlations Li, behave algebraically when the 


(e(y-P) 1)2’ 
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chemical potential is in the bulk band whereas they have 
combinations of algebraic and exponential dependence 
when p. is in the bulk gap. Meanwhile, the figure of merit 
is larger when the chemical potential is away from the 
band edge. 

Since the actual thermoelectric transport coefficients 
involve both the edge and bulk correlation functions, we 
use = L® + L^, and compute the transport properties. 
It has already been mentioned in the earlier section that 
the thermoelectric transports in the topological systems 
are due to the competition between the edge and the 
bulk contributions. As has already been discussed 
above, the relative magnitudes of both the edge and 
the bulk transports for different chemical potentials 
are different, the figure of merit ZT is also different. 
When the chemical potential is in the bulk band, ZT 
from the edge becomes larger eventhogh the number of 
edge carriers are exponentially small. This is due to the 
fact that that the edge transport is overridden by the 
bulk contributions and a high ZT from the edge never 
appears. In a similar way, ZT from the bulk becomes 
larger when the chemical potential is in the bulk gap. 
In this case, there are very few bulk carriers whereas 
there are some edge carriers for which the edge carriers 
dominate and hence ZT is suppressed. This is a clear 
indication that in topological system the edge and the 
bulk transport compete with each other suppressing the 
total ZT. 

In order to compute the transport properties, we con¬ 
sider the following parameters: the hopping integral is 
taken to be zero and the next nearest neighbour hop¬ 
ping integral as t' = —0.3z. The SC and SDW order 
parameters respectively are considered to be Aj = Ag = 
0.0050, Mq = 0.0750 where the parameter 0 taken to 
be 0 = O.SOeC. The computation is performed at a 
temperature T = 20Ar. The relaxation time is taken 
to be Te = I0“^^sec. The phonon thermal conductivity 
is taken to be constant which is kl = 0.5Wm~^K~^. 
The effective system size £ is assumed to be I/xm and 
the cross-sectional area s is taken as lOnm x O.lnm. All 
these parameters are taken in an adhoc basis and the 
thermoelectric coefficients are computed. The results are 
shown in Fig. 2. and Fig. 3. It is obvious from Fig. 2 
and Fig. 3 that the edge contribution to ct, Kg and ZT is 


much higher as compared to the bulk. On the contrary, 
the Seebeck coefficient has comparable contribution from 
both, a positive contribution from edge whereas a neg¬ 
ative contribution from bulk. Thus, there is a partial 
cancellation due to opposite charge carriers in the edge 
and the bulk. This results a peak structure in the Seebeck 
coefficient near the band edge which makes the figure of 
merit large. 

Due to the edge-bulk competition, there is a possibility 
that ZT could be maximum when the chemical potential 
is near the band edge. In such a case, the bulk conduction 
is dominant at high temperature. On lowering temper¬ 
ature, the bulk-edge cross-over takes place and there is 
a possibility that ZT might be higher as shown in Fig.2 
and Fig.3. The edge current in these system provides bal¬ 
listic transport but such transport crucially depends on 
the inelastic scattering length £ of the edge states because 
they loose their coherence due to inelastic scattering. At 
high temperature, £ is short and the bulk transport is 
dominant but at low temperature, £ is long and the edge 
transport is dominant. This might be the reason why an 
enhanced figure of merit results at low temperature. 

In conclusion, we summarize the main findings 
of the present work. We study the thermoelectric 
transport in a two-dimensional topological system 
which has coexistence of superconductivity(SC) and 
spin-density wave(SDW). The SC is presumed to be 
of dx 2 -y 2 + {px + ipy) type whereas the SDW order 
parameter is of BCS symmetry. The Hamiltonian 
having such a structure is shown to have topological 
coexistence phases in addition to the conventional one. 
The thermoelectric transport in such topological system 
is discussed with respect to two distinct contributions: 
(i) the surface and (ii) the bulk. The competetion 
among the surface versus the bulk transport is analyzed 
in different parameter regimes and the possibility of en¬ 
hancement of figure of merit is discussed in these systems. 
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£r(l/(/jflm)) 






FIG. 2. Computation of longitudinal electrical conductivity, 
Thermopower, electronic thermal conductivity and figure of 
merit as a function of the chemical potential (/r) at T = 20K 
for £ = l^m. Black curve shows the bulk, dashed curve shows 
the edge and red curve shows the total contribution to the 
above quantities. For illustration, here, we have chosen Aj = 
Ag = 0.005«, Mo = 0.075a, t' = -0.3a, kl = 0.5Wm-^K-\ 
c = O.lnm (a = 0.30 eV). 




















7 


£r(l/(/jflm)) 






FIG. 3. Computation of longitudinal electrical conductiv¬ 
ity, thermopower, electronic thermal conductivity and fig¬ 
ure of merit as a function of the chemical potential(/r) at 
T = 2QK for £ = 0.8^m. Black curve shows the bulk, dashed 
curve shows the edge and red curve shows the total contri¬ 
bution to the above quantities. For illustration, here, we 
have chosen Aj = Aq = O.OOSz, Mo = 0.0752, t' = —0.32, 
kl = Q.hWm~^K~^, c = O.lnm (2 = 0.30 eV). 




















